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We study lepton number violating decays of charged K, D, D s , B and B c mesons of the form 

M + — > M' £ + £ + , induced by the existence of Majorana neutrinos. These processes provide 

information complementary to neutrinoless double nuclear beta decays, and are sensitive to 

£"S ' neutrino masses and lepton mixing. We explore neutrino mass ranges m^ from below 1 eV to 

\Q ' several hundred GeV. We find that in many cases the branching ratios are prohibitively small, 

however in the intermediate range m n < m^ < tub , in specific channels and for specific 

neutrino masses, the branching ratios can be at the reach of high luminosity experiments 

like those at the LHC-6 and future Super flavor-factories, and can provide bounds on the 

C^ ' lepton mixing parameters. 
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One of the outstanding issues in neutrino physics today is to clarify the Dirac or Majorana 



character of neutrino masses. The discovery of neutrino oscillations indicates that neutrinos are 

o: 

massive particles with masses likely to be much smaller than those of charged fermions [l|. This fact 
provides an important clue on the existence of a more fundamental physics underlying the standard 
S^ model (SM) of particle physics, because neutrinos are naturally massless in the SM. Although the 

experimental results on neutrino oscillations can determine the neutrino mixing parameters and 
their squared mass differences, the absolute magnitudes of the masses as well as their origin remain 
unknown and constitute fundamental open questions in neutrino physics. Many experiments have 
been set to search for the absolute magnitude of neutrino masses. Direct methods to determine 
the mass of the electron neutrino use the endpoint of the electron spectrum in beta decays. The 
most sensitive of these experiments uses Tritium [2J], setting the present upper bound m Ue < 2 eV 
[3], and the next experiment is expected to reach a sensitivity of 0.2 eV [4]. Other experiments do 
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direct searches for muon and tau neutrino masses, setting the upper bounds m v < 190 keV and 
m VT < 18.2 MeV respectively, at 90% C.L. [3(. To date, the most stringent bound on the sum of 
all light neutrino masses is obtained from cosmological observations, given by Y2i m vi < 0-1? eV 
at 95% C.L., a figure which is, to a certain extent, model dependent [5]. 

If neutrinos are Dirac particles, they must have right-handed electroweak singlet components 
in addition to the known left-handed modes; in such case lepton number remains as a conserved 
quantity. Alternatively, if neutrinos are Major ana particles, then a neutrino is indistinguishable 
from its antiparticle and lepton number would be violated by two units (AL = 2) in some processes 
that involve neutrinos. The experimental results to date are unable to distinguish between these 
two alternatives. 

There have been several attempts to determine a Majorana nature of neutrinos by studying 
AL = 2 processes. The most prominent of these processes are neutrinoless nuclear double beta 
decays (Ov/3/3), which have been regarded as the most sensitive way to look for lepton number 
violation (LNV) |6(]. The observation of Oz^/3/3 would indeed be very important not only because 
it would establish the existence of LNV - implying that neutrinos are Majorana particles, but 
also because they would provide a scale for the absolute magnitude of light neutrino masses, 
complementary to the direct searches mentioned above: these nuclear processes are proportional 
to the square of the effective neutrino mass m ee = \ Yli=i ^e« m *l> w hh m, and U e i being the 
individual neutrino masses and the vi — e mixing matrix elements, respectively [3|- However, it 
has long been recognized that, even though the experiments are very sensitive, the extraction of 
the neutrino mass scale and the Majorana nature of neutrinos from nuclear Ov/3/3 is a difficult 
task, because reliable information on neutrino properties can be inferred only if the nuclear matrix 



elements for 0^/3/3 are calculated correctly. The calculation of the nuclear matrix e 
usually performed within either the quasi-particle random phase approximation 



ements for 0^/3/3, 
or the nuclear 



shell model [9[] or their variants, is known to be a complex task, sometimes with large differences 
among the different approaches (lOJ. Even in the most refined treatments, the estimates of the 



nuclear matrix elements remain affected by various large uncertainties 
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Another avenue to detect the Majorana nature of neutrinos is to study AL = 2 processes in 



rare meson decays 
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14]. In this paper we study AL = 2 decays of heavy charged mesons whose 
signals could be captured at high intensity experiments such as LHC-6 and future Super L?-factories 
as well as advanced LT-factories. The AL = 2 processes we treat in this paper are rare neutrinoless 
decays of heavy charged mesons into a lighter meson and two charged leptons of the same sign [13(] . 
These processes, just like neutrinoless nuclear double beta decays, can occur only via Majorana 



neutrino exchange, and thus their experimental observation could establish the Majorana character 
of the neutrinos and the absolute scale of neutrino masses in much the same way as in nuclear Qu(3(3 
decays, but there are some essential differences. Prom a theoretical viewpoint, the uncertainties 
in meson decays are much easier to handle than in nuclear 0vf3f3 decays. However, from the 
experimental viewpoint, the AL = 2 meson decay rates in the case of standard neutrinos (m v < 2 
eV) are prohibitively small for any experiment, while Oz^/3/3 decays are more realistic options, due 
to their macroscopically large samples of decaying nuclei. In contrast, for heavier, non-standard, 
neutrinos, the meson decay rates are good alternatives to search for, as they can be within reach 
of future experiments. 

In this study it is important to distinguish between standard and sterile neutrinos. From direct 
searches we know the standard electron neutrino mass is below 2 eV [5|], and neutrino oscillation 
experiments tell us that all three neutrino masses differ from one another by much less than that 
value [15l |. Therefore all neutrinos with masses above 2 eV are assumed to be non-standard. Since 
our work is mainly relevant for neutrinos above this bound, in what follows we will denote them 
generically by the letter N, instead of v. 

An important motivation to search for sterile (non-standard) neutrinos with masses of the 

order of 1 MeV is that their existence has nontrivial observable consequences for cosmology and 

astrophysics. They are presumed to participate in big-bang nucleosynthesis, supernovae explosions, 

arge scale structure formation and, in general, to be a component of the dark matter in the universe 



161 ] . Thus, sterile neutrino masses and their mixing with the standard neutrinos must be subjected 
to cosmological and astrophysical bounds [17]. There are also some laboratory bounds coming from 
the fact that sterile neutrinos contribute via mixing with the standard neutrinos to various processes 
which are forbidden in the SM. Those bounds turn out to be much weaker than the cosmological 



and astrophysical bounds, but useful in cases where the latter become inapplicable [18]. 

We have separated the analysis into three different cases, depending on the relevant neutrino 
mass range. If the exchanged neutrino is much lighter than the energy scale in the process, the 
amplitude of the decay rate is proportional to the square of an effective electron-neutrino mass, 
TO ee = I ^2n UeN m N\ 2 , which is anticipated to be of the order ~ 1 eV 2 or less from current neutrino 



data and cosmological observations such as WMAP 19], if only standard neutrinos are involved. 
Instead, if the exchanged neutrino is much heavier than the decaying meson, the decay rate is 
proportional to \UN£ 1 UNi 2 / m N\ 2 , where m^ and Uni are the heavy neutrino mass and its mixing 
with the standard leptons, respectively. In general, in this case U^i is small and tun is large, so 
the factor constitutes a severe suppression to the decay rate. Finally, for the case of Majorana 



neutrinos with intermediate masses between that of the initial and the final meson, the decay rate 
is dominated by a resonantly enhanced s-channel amplitude 
neutrino goes on its mass shell. 
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where the intermediate 



In Section II, we describe the approximation methods for the calculations of rare heavy meson 
decays of the form M + — > M'^lfl^ (where M and M' are pseudoscalar mesons). Here we are 
interested on K + , D+, D+, B + and B+ decays into tt~£ + £ + , K~£ + £+, D~£ + £ + , D~£ + £ + and 
B~£ + £ + , (where £ = e, fj, or r), therefore, we will denote the initial and final mesons generically 
by M + and M'~ , respectively. We separate the analysis for the three cases of light neutrinos 
(mjv < TnM')i intermediate neutrinos {tum 1 < w^N < ^Af), an d heavy neutrinos {m^ > tum)- We 
include the results and discussions in each subsection. In Section III we summarize the results and 
state our conclusions. 

II. CALCULATIONS OF M+ -> M'-l+l+ 



We now describe our approximation methods for the calculations of rare heavy meson decays 
of the form M + — > M'~lfl2 (where M and M' are pseudoscalar mesons) in all three neutrino 
mass ranges described above. At the quark level, the decay occurs via two types of amplitudes, 
shown in Fig. [TJ We find that in the case of light neutrinos {niN < Tn n ), the amplitude on the 
left in Fig. [T] ("t-type" diagram) dominates due to long distance contributions, and the decay rate 
becomes proportional to m 7 M x m 2 N . For this reason, only the decays of the heavier B mesons are 
of any importance in this case. In contrast, for intermediate neutrino masses {jum 1 < mN < v ^m), 
the diagram on the right in Fig. Q] ( "s-type" diagram) dominates when the neutrino propagator 
becomes resonant on its mass shell, in which case the decay rate turns out to be less dependent 
of the neutrino mass, but very sensitive to the mixing elements. Finally, for heavy neutrinos 



(m-Tv > thm), both amplitudes in Fig. Q]are comparable and the decay rate is oc 1/m 
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FIG. 1: The i-type and s-type weak amplitudes at the quark level that enter in the process M^ 
(plus the same diagrams with leptons exchanged if they are identical) . 
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TABLE I: Values of input parameters used in our calculations. They correspond to the central values given 
in Ref. [3(, except for fg and fg t which are taken from Ref. 2l|, and V cs which is calculated by imposing 
the unitarity constraint on the CKM matrix. 



Parameter 


Value 


Parameter 


Value 


h 


130.4 [MeV] 


V ud 


0.9742 


Ik 


155.5 [MeV] 


v us 


0.2255 


fLH- 


205.8 [MeV] 


v ub 


0.0039 


fn a 


273. [MeV] 


v cd 


-0.230 


Ib 


196. [MeV] 


Vcs 


0.950 


/*. 


322. [MeV] 


v cb 


0.041 



In TableUwe list the numerical values of the input parameters we use in our numerical estimates. 



A. The case of light neutrinos (tojv < rn w ) 



We find that a neutrinoless decay like B + — > D~t + £ + with light Majorana neutrinos in the 
intermediate state is dominated at the meson level by the amplitude shown in Fig. [21 when the 
intermediate state goes on mass shell. This amplitude originates at the quark level from the t-type 
weak amplitude shown in Fig. [H We find the s-type amplitude shown in Fig. [T]to be subdominant, 
or at most comparable with the former. In this sense, our treatment differs from that of A. Ali et 
al. [131 ]. where the s-type amplitude is assumed to dominate 22( | . However, since the rate in any 
case turns out to be too small for any foreseeable experiment, we will just do an order-of-magnitude 
estimate for it, calculating the absorptive part and assuming that the dispersive part is not much 
larger. The absorptive part of the amplitude is calculated by setting the intermediate particles on 
their mass shell and then integrating over their phase space: 

(1) 



-Mabs(^ ->■ D rr) = / dps DN A B -n) N i A DN ^ m 

where Ab^dn£ and Adn^-di are the tree- level amplitudes for the respective sub-processes, and 
dps DN is the Lorentz-invariant phase space of the intermediate D-N pair, which in the rest frame 
of the pair is dps DN = ^ s (l/Wir 2 )(\p N \/m£)£)dQ,N- Here ^ s is the sum over the neutrino spins, 
Piv is the 3-momentum of the neutrino in the D-N rest frame, and mot is the invariant mass of 
the pair. In turn, the amplitudes of the weak sub-processes are: 
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FIG. 2: The main diagram in an effective meson theory for M + — > M'~£ + £ + (plus diagram with leptons 
exchanged if they are identical), mediated by Majorana neutrinos, when the neutrino is much lighter than 
the final meson. The amplitude is estimated considering the intermediate state on its mass shell. 

A B +-+D°m = % Vcb U Nt (D°(p')\J»(0)\B + (p)) u N ( PN )^(l - l5 )v e (h) (2) 

A D o N ^ D - e = ^V ud U m (D~(p')\J»(0)\D°(p))v N (p N ) lfl (l- l5 )v e (l 2 ), 

where V is the Cabbibo-Kobayashi-Maskawa (CKM) matrix for quark mixing, and U the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix for lepton mixing. The hadronic matrix el- 
ements can be parameterized in terms of phenomenological form factors F + {q 2 ) and F~(q 2 ) as 

(D (pi)\J»(0)\B+(p))=F+ D ( q 2 )(p + pr + F BD ( q 2 )(p- P r, (3) 

and similarly for (D~(p')\ J^ (0)\D° (p)) , where q is the corresponding 4-momentum transfer. In our 
crude estimate, we will neglect the F~ form factors and assume the F + to be constants of order 
unity over the kinematical range. 

Now, in the product Ab^-dni x Adat->D£, after summing over intermediate spin states, the two 
lepton lines can be combined into a single one by using in Adn->D£ the identity ^7^(1 — 75)^ = 
U|7^(l + 75)uat. Here Ug is a -u-spinor for the charged antilepton, and the neutrino is assumed to 
satisfy the Majorana condition u^ = Xn un (where Xn is a phase). The result is then: 

M abs = ^V cb V ud U 2 m X N I ^Mf+ dF + d (4) 

ui(h)(^D+ $D°)( l + 7s)(^iV + m N )(p' D o+ #b)(1 - 75)^(^1), 

This angular integral is quite simple, because in the D-N frame the energy of every particle in the 
process is fixed. The subsequent steps to obtain the decay rate are straightforward and described 
in the Appendix. The expression for the rate is thus the integral [see Eq. (|A8|) ]: 



(tub— mi) 

g f i T/ T/ | 2f+2 F+2 \Um^LNi ] f dmpj |p r ' 2 

(16^)2 I V cbVud\ t BD t DD -g- J — — ^ 

(m D +m e ) 



ns + -► B-tt-) = t^IVMTf+Zf+I >"Z7 ' / ^ ^-\h\ II2I x n, (5) 



where |pat|, |1i| and | X^ | are the 3- momenta of the neutrino and leptons (given in the Appendix) and 
1Z is a quantity of dimension m 6 shown in Eq. (|A6P . The integral can be easily done numerically, 
which we do considering a D meson in the intermediate and final states (b — > c transition), or 
alternatively a pion (b — > u transition). 

Notice that by assuming the form factors to be constant unity we are overestimating the process, 
while by neglecting the F~ form factors and the dispersive part of the amplitude we may be inducing 
an uncertainty of an order of magnitude. Within our approximations, in both cases the results for 
the branching ratios are extremely small: 

Br (B + -> D-£+£ + ) ~ 1.2 X 1(T 31 (^f^f) , (6) 

Br {B + -»• 7t~l + e + ) ~ 2.3 x 1(T 33 C^f^f) , (7) 

where we used the values of the CKM elements shown in Table HJ and also Fb = 4.0 x 10 -13 GeV. 
We can compare these results with those of A. Ali et al. [13], who considered the s-type diagram 
only. In our notation, their result for Br(B + — > 7r _ e + e + ) becomes (0.3—1.8) x 10 -35 (Cjy»mjv / eV) 2 , 
which is two orders of magnitude smaller than Eq. (J7|) . 

Nevertheless, we expect our results to be just rough estimates within one or two orders of 
magnitude, as we have taken the form factors Fg D ~ Fnn ^° ^ e un ityj and we have neglected the 
form factors F^ D and F^ D altogether. In general, the form factors F + are expected to be unity at 
most at the kinematical end point where the two meson wave functions could overlap completely 
(provided they have the same shape), but it should be smaller for all other q 2 values. 

Taking for F + an average value of e.g. 0.3 instead of unity, our calculated rates get reduced by 
a factor (F + ) 4 ~ 10~ 2 , reducing Eq. ([7]) to a value comparable with the result of Ali et al. 

Accordingly, in the case of light neutrinos, our crude estimate cannot clearly show the dominance 
of the t-type diagram. However, it does show at least that a calculation based purely on the s- 
type diagram may be an underestimation [221 ] . It also shows that this potential underestimation is 
hardly more than two orders of magnitude, keeping these branching ratios still beyond the reach 
of foreseen experiments, as concluded in Ref. [131 ]. 

To estimate the actual range of these branching ratios we would need to have estimates of the 
neutrino masses and mixings as well. Using the standard parametrization of the PMNS neutrino 
mixing matrix multiplied by a 3 x 3 Major ana phase matrix, the term U^/tun can be explicitly 
written in terms of three light neutrino masses, three neutrino mixing angles, two Majorana phases 
and one Dirac phase. Since the sign of AmJ^ is not determined from the existing data, there 



are two possible neutrino mass hierarchies, one called normal (771,3 > ^1,2) and the other inverted 
(771.3 < "^1,2)- The size of the term Uf^fn^ in general depends on the mass hierarchy. If we consider 
standard neutrinos, we know that m u < 2 eV, and we can roughly use U U £ ~ 0(1) for either £ = e 
or fi, in consistency with oscillation experiments. We then get branching ratios smaller than 10~ 31 
and 10 -33 , respectively, values which are prohibitively small for any foreseen experiment. On the 
other hand, if we consider heavier neutrinos (but still lighter than m n ), i.e. m^ ~ 100 MeV, the 
results could be more promising, but in those cases we should use the mixings of standard with 
extra neutrinos, which are suppressed: Ufj e ,Uff < 0.002, [23(] so the resulting branching ratios 
have the upper bounds 10 -21 and 10~ 23 , respectively, which are still prohibitively small. 

As a final remark, we want to comment on the assumptions involved in this calculation. First, 
the fundamental process at the quark level (see Fig. 1) with two electroweak vertices has been 
modeled as a process with hadrons and leptons, where a single long distance contribution (an 
intermediate state with a meson and a neutrino on shell) is supposed to dominate; we have thus 
neglected other possible intermediate hadronic states (e.g. excitations of the intermediate meson 
and multimeson states) as well as a short distance contribution where both weak vertices coalesce 
into a single one [24] . We have assumed the dominance of the single D-N intermediate channel as 
it goes on its mass shell. Another assumption was to consider the absorptive part as representative 
of the full amplitude; since we are only after an order of magnitude estimate, this is likely to be 
a good assumption, again due to the resonant character of the intermediate state as it goes on its 
mass shell. Within the hadronic approximation for the weak currents, we took into account just 
one of the form factors of each hadronic current, and assumed it to be constant (unity) within the 
whole dynamic range. In principle one can expect the form factor to be unity at most, as explained 
before; taking the q 2 dependence into account one should then obtain a lower value for the rate, 
but as we have seen, it is unlikely for this effect to change the result by more than two orders of 
magnitude. These approximations are therefore consistent with the level of precision we seek. 

B. The case of intermediate mass neutrinos {m n < m jy < msj) 

In contrast to the previous case, the process M + — > M' £ + £ + in the case of Majorana neutrinos 
with masses in the intermediate range txim' < rriN < thm is dominated by the s-type amplitude of 
Fig. [H corresponding at the meson level to Fig. [3J as the neutrino in the intermediate s-channel 
goes into its mass shell. As stated in the Introduction, Majorana neutrinos with such masses must 
be sterile and should originate from new physics beyond the SM. 
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M + M" 

FIG. 3: The dominating diagram (plus diagram with leptons exchanged if they are identical) in an effective 
meson theory forM + — > M' £ + £ + , mediated by Majorana neutrinos with mass in the range between tum' 
and rriM- 

Since there are two identical leptons in the final state, one must also consider the diagram 
with crossed leptons and then integrate over half the phase space. However, for the case of the 
intermediate neutrino on mass shell, the result is the same as using a single diagram, as if the 
leptons were distinct, as shown in Fig. El The effective amplitude then is: 

•M = -TT U N£ V* Q V* q f M fM' 7-3 2 \ , ■ F~ ^ 

2 « (p* N - m z N ) + im N T N 

where Uni and Vq iq are the PMNS lepton mixing and CKM quark mixing elements, respectively, 
/a/j Im' are the meson decay constants, and we define M as the reduced matrix element that 
contains all the spinor structure of the amplitude: 

M = \ N %(Zi)j4r(l +75) (&N + m N )rf M >(l -75)^2) (9) 

where the notation is the same as in Eq. Q. 

The decay rate we seek is then given by T{M + — > M' £ + £ + ) = (l/2mjw) fdps 3 \A4\ 2 , where 
<ips 3 is the final 3-particle phase space. The calculation of the squared matrix element and the 
integration over the final phase space are shown in Appendix 2, resulting in the following expression 
[see Eq. flB3J) ]: 



vy f/„2 , 2n 2 I 2 2\2~l (/ 2 2\2 / 2 , 2\ 2 ~| 

x {(m N + m t )m M - (m N - m t ) ){{rn N -m l ) - (m N + m e )m M ,} , 



where | li | and 1 12 1 are the 3- momenta of the first electron in the M meson rest frame and of the 
second electron in the neutrino rest frame, respectively. 

Before we can use this expression, we also need a theoretical expression for Ttv, the total decay 
width of the intermediate Majorana neutrino, in terms of the same neutrino parameters we have 
just used. The total width Ttv can be estimated by comparing the decay modes of N with those 
of the t~ lepton, where T T oc m^.. Both N and r~ decay via the same type of diagrams and 



10 

couplings, but there are a few differences: (a) N has a different mass (thus Tjv oc m 5 N ); (b) Tjy 
has an additional factor of two due to the Majorana character of N (unlike t~ which is a Dirac 
particle), because it decays with equal probability into both [£'~ +rest + ) and (£' + + rest ); (c) IV 
has an additional mixing factor | t/jv^' 1 2 - Therefore: 

5 



m 



^ N ^2^\Um>\ 2 {^) xT T , (11) 



I' 

This expression for T^r is a good approximation when m^ is near 2 GeV; in this case the decay 
channels of N are those of r, where the virtual W boson produces e~V e ,n~Vu, and du (the last 
channel is actually a set of three, due to color). However, for tun > 2 GeV, the additional channels 
t~V t and sc open, increasing the expression in Eq. (jlip by up to a factor ~ 1.5, including phase 
space suppression due to the masses of the products. Consequently, using Eq. ([11]) in Eq. (flQj) may 
oversetimate the rates by at most ~ 30%. We will thus use Eq. (Ilip in the estimation of the LNV 
rates, but keeping in mind that a correction in Ttv should be included in a more refined study. 
Accordingly, and if we neglect the charged lepton mass, Eq. (|10p turns into: 

,5 / ■>■ s 2 



r ( M - mw) - lauyuv^f J^ . mM - m * d - =*Hi - z 



lit 

■2 



12fk2 -*■•/«,«-,. «■«.«, & |£, w ,|2 2 r x V - m5r ; V - ^ 

(12) 

Here we will use m r = 1.77 GeV and T T = 2.3 • 10 12 GeV [3J. Eq. (112j) is valid for rajy in the range 



m^' < tun < mMi it vanishes at the two endpoints of this range, and reaches its maximum at 
m N = yjrriM -tum', where (1 - m 2 M ,/m 2 N ) 2 (l - m 2 N /m 2 M ) 2 -)■ (1 - m M '/m M Y- 

Consequently, these suppressed non-standard decays can impose more or less stringent bounds 
on the mixing elements between the standard leptons and extra neutrinos, |[/jvd> depending on 
the Majorana neutrino mass. In particular, the non-observation of these processes defines tun- 
dependent upper bounds for the corresponding |C/jv^|- 

In Figs. [3H3 w e show the branching ratios for the decays K + — > tt^£ + £' + , D + — > M'~£ + £ + , 
D+ -> M'~£ + l + , B + ->• D~£ + £ + and B+ -)■ M''£ + £ + as functions of m N , where the bounds on 
the mixings \Uni\ can be deduced also as functions of m^- 

Let us consider the decay B + — > D~e + e + as an example. Here we must use V q Q — > V u b and 
V qiq2 — > V c d as inputs, as well as fs and fjj (see Tabled]) and Tb = 4.0 • 10~ 13 GeV [3j. The 
branching ratio for this process as a function of mjy is shown in Fig. EKa), lower dashed line, and 
reaches a maximum: 

Br max ( J B+ -> D~e + e + ) = 3 • 10~ 7 x ' ^ e| at m N ~ 3 GeV. (13) 

2^/ l"w| 
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FIG. 4: Branching ratios for K + — > it £ + £ + (£ = e, fj) as functions of the exchanged neutrino mass mjv in 
the range m^ < mjv < rriK, with the lepton mixing factor, |t^Ve| 4 /X^' I t^JV£' | 2 , divided out. 
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FIG. 5: Branching ratios for (a) D + decays and (b) D+ decays, as functions of the neutrino mass thn, with 
the lepton mixing factor divided out as in Fig. SJ The full lines correspond to I = e and the dashed lines to 

This expression just gives the maximal possible value of this branching ratio, which occurs only if 
rriN happens to be near 3 GeV, but for other values of win, it could be much smaller, as shown in 
Fig. Ha). 

Analogous to Eq. (I13D , the maximal branching ratio of any of the other decays has the form: 

\U m \ A 



Br max (M + -► M' 



-)=Cx 



(14) 



J2e' \ U N£'\ 2 ' 

Table HI1 shows the coefficient C appearing in Eq. (|14p . for the different branching ratios, and the 
value of the corresponding neutrino mass m^ at which the maximal branching ratio is reached. 

Accordingly, an experimental upper bound on the branching ratio for M + — > M'~l + £ + imposes 
an upper bound on the leptonic mixings |£7jv^|j bound that strongly depends on the neutrino mass 
ttitv , and which is most stringent if m^ ~ \Zttim ■ tum'i where the branching ratio is maximal. For 
win away from that value, the upper bounds imposed on the mixings become much less stringent. 

From the C values in Table lU one can read the potential of different processes to set upper bounds 
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FIG. 6: Branching ratios for B + — > M'~£ + £ + as functions of the neutrino mass mjv, with the lepton mixing 
factor divided out as in Fig. |H The produced pseudoscalars are M' = n, K, D, D s . (a) The case of leptons 
with negligible mass (£ — e,fi); (b) the case £ = r (here M' = D,D S are kinematically forbidden). 
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FIG. 7: Branching ratios for B c — > M'~£ + £ + as functions of the neutrino mass rajq, with the lepton mixing 
factor divided out as in Fig. |4j The produced pseudoscalars are M' = n, K, D, D s . (a) The case of leptons 
with negligible mass (£ — e, \x)\ (b) the case £ = t. 

on the lepton mixing elements | C/jv^ I » f° r different neutrino masses m^. For a given experimental 
upper bound of a branching ratio, the larger the C coefficient, the more stringent the upper bound 
that can be imposed on |E7jv^|, provided the neutrino mass is near the indicated value where the 
theoretical branching ratio is maximal. 

From Eq. (I14p it is clear that the bounds on the mixings imposed from these decays appear in 
the combination 

\U m \ A 



e, jj, or r, 



(15) 



|t/iVe| 2 + |^| 2 + |^iVr| 2 ' 
not just |£/jvi|- Only if \Un(\ is much larger than the other mixings, then this expression reduces to 
\Uni\ 2 ■ Otherwise, one must use the bounds on Br(M — > M'~ £ + £ + ) for a given meson pair M and 
M', but for all lepton flavors £ = e, /jl, r, in order to disentangle the bounds for each of the mixings 
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TABLE II: The coefficients C appearing in Eq. (|14p for the maximal branching ratio, and the neutrino mass 
tun at which the maximum is reached, for various decays M + — > M'~£ + £ + , where mg can be neglected. In 
the last column, the expected upper bound on the branching ratios, provided |£//vf| 2 ~ 10 -6 or IO -7 , for 
niN ~ 0.1 GeV or ~ 1 GeV, respectively. 



decay 


C 




mjv at maximum 


Br < 


K+ -> tt-£+£+ 


2.8 




0.26 GeV 


2.8 


io- 6 


D+ -> 7T-£+£+ 


4.5 


10- 3 


0.51 GeV 


4.5 


io-i° 


D+ -> K-£+£+ 


1.4 


io- 4 


0.96 GeV 


1.4 


IO-" 


D+ -> ir-£+£+ 


6.9 


io- 2 


0.53 GeV 


6.9 


io- 9 


D+ -> K-£+£+ 


2.2 


IO- 3 


0.99 GeV 


2.2 


io-i° 


D+ -> D-£+£+ 


8.5 


io- 8 


1.92 GeV 


8.5 


io- 15 


B+ -> 7r-£+£+ 


6.3 


io- 6 


0.86 GeV 


6.3 


IO-" 


S+ -> A'-^+^+ 


3.6 


io- 7 


1.61 GeV 


3.6 


IO-" 


S+ -> D-£+£+ 


1.7 


io- 7 


3.14 GeV 


1.7 


10-" 


B+ -> D s 7£+£+ 


4.5 


io- 6 


3.23 GeV 


4.5 


10-13 


B+ -> 7r-£+£+ 


6.4 


io- 4 


0.94 GeV 


6.4 


IO-" 


B+ -> jr-^<-£+ 


3.9 


io- 5 


1.76 GeV 


3.9 


10-1 2 


B+ -> D-£+£+ 


2.4 


io- 5 


3.43 GeV 


2.4 


10-12 


B+ -> ^ s 7£+£+ 


6.5 


io- 4 


3.52 GeV 


6.5 


IO-" 


B+ -4 B-£+^+ 


1.6 


IO-" 


5.76 GeV 


1.6 


10-1 8 



|C7jv^|. Moreover, these bounds will depend on m^, since the relation between the branching ratios 
and the mixings depend on tun, as it was already mentioned and shown in Figs. SHZl 

On the other hand, to explore the prospects of experimentally observing any of these processes, 
one needs at least an estimate of the \Uni\ elements. Present upper bounds on the heavy-to-light 
neutrino mixing (L^l 2 for £ = e,/x, vary considerably with the neutrino mass, but are typically in 
the range \Um\ 2 < 10" 4 , 10" 6 , 1(T 7 , for m N ~ 10 MeV, 100 MeV, 1 GeV, respectively (pp. 546- 
548 in Ref. 3j). We have then listed in the last column of Table ITT1 the expected upper bound on 
the branching ratios, provided the mixing elements have the values just mentioned. 



C. The case of heavy neutrinos (tun > rris c ) 



If neutrinos happen to be much heavier than the decaying meson, then in general both diagrams 
in Fig. 1 contribute with more or less the same strength, and reduce at the meson level to a single 
point-like interaction diagram as shown in Fig. [HJ The vertex in Fig. [8] represents the double weak 
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FIG. 8: The diagram in an effective meson theory for M + — > M' e + e + , when the neutrino mass is much 
larger than that of the decaying meson, which is resulted from the four amplitudes of Fig. 1. 

interaction shown in Fig. 1, where the neutrino line as well as all other internal lines have been 
reduced to a point. At the meson level, the specific tensor structure of this four-particle vertex 
cannot be selected among all the general possibilities, so we start from the fundamental quark and 
lepton interactions as shown in Fig. 1 and exhibit the approximations involved to get to the leading 
term at the meson level. These details are presented in Appendix C. Our model of the dynamics 
in this case is equivalent to that of Ali et al. [13[ • In summary, if we can approximate the hadronic 
tensor by the product of two currents, factorized by a vacuum insertion, the squared amplitude 
is then given in terms of the mesons' decay constants and the kinematics of mesons and leptons 
separate into independent factors (see Eq. IC5|) : 



\M\* ~ fttfoipM-pM'Yitl-h). 



The decay rate then becomes (see Eq. IC6h : 



r(M + 



-»• M'-£+£+) 

(m M -m M ,) 2 



128vr 3 



7-7*2 



m N 

2 



dm\ t ^\l^ '"" 



Ami 



2 ' 2 
m M m M 



V*nV* + 

iQ 9192 ~ 

A x / 2 (1 



9l<2 992 



N r 



/at/at m M 



2 2 

m| mi 

' 2~' 2~ 

m u rn u 



) 1 + 



ni 



M' 



ml. 



in 



M 



m 



M 



(16) 



\m, 



2mi 



where the function \(x,y,z ) is defined in Eq. (|A2j) . This expression exactly coincides with the 
expression obtained in Ref. 13[] for the heavy neutrino cases. 

The integral above can be easily done numerically. In order to do the phenomenology, we 
set a fiducial value for the neutrino mass rriN = 100 GeV, and a corresponding mixing element 
I ^Jnh 1 2 = 10 -2 and express the branching fraction of this decay in terms of a dimensionless quantity 



B, whose value, according to Eq. (|16p . is determined by the masses of the external particles: 



For the case B + — )■ D £ + £ + we must use fu = Ib+ an d fw = Id- ■ Using the values shown 
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TABLE III: Branching ratio coefficients B appearing in Eq. (fTTj) . for various decays M + — > M'~£ + £ + , if the 
process is dominated by heavy neutrinos (tun 3> rnu)- & values correspond to branching ratios if mjv = 100 
GeV and |£/jvf| 2 = IO -2 . All three lepton flavors are considered (£ = e,fi,r). Entries are empty for decays 
that are kinematically forbidden. 



decay 


B(£ = e) 


B(* = /i) 


B(£ = t) 


K+ -> ir-£+e+ 


8.47 


io- 24 


2.44 • 10- 24 


- 


D+ -» 7r-£+£+ 


1.90 


10 -23 


1.78 -IO- 23 


- 


D+ -> K-£+£+ 


1.58 


10 -23 


1.47 -IO- 23 


- 


D+ -» 7r-^+^+ 


2.14 


1Q -22 


2.02 • IO" 22 


- 


L>+ -> K-£+£+ 


2.46 


10 -23 


2.30 • IO- 23 


- 


D+ -> D-£+£+ 


6.99 


10 -28 


- 


- 


B+ -> 7r-£+£+ 


1.13 


1Q -23 


1.12 -IO- 23 


7.42- IO" 25 


S+ -> A'-^+^+ 


8.44 


10 -25 


8.37 -IO- 25 


5.01 • IO" 26 


S+ -> D-£+£+ 


1.02 


1Q -22 


1.01 • io- 22 


- 


B+ -> Dj^+£+ 


5.02 


10 -23 


4.96 • IO- 23 


- 


B+ -> 7r-£+£+ 


1.76 


io- 21 


1.75 -IO- 21 


3.04- IO" 22 


s+ -> i?-^+£+ 


1.73 


10 -22 


1.72 • IO" 22 


2.89- IO" 23 


S+ -> D-£+£+ 


3.20 


1Q -22 


3.17 -IO- 22 


2.14- IO" 23 


B+ -> D-£+£+ 


9.17 


io- 21 


9.10 -IO- 21 


5.17- IO" 22 


B+ -> B-£+£+ 


3.31 


io- 28 


2.96 -IO" 28 


- 



in Table U as well as Y B + = 4.0 • 10 
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the result is: 



Bt{B + -*■ D~£ + £ + ) rj 1.1 • IO" 22 



100 GeV 

tun 



rm 2 - 



(18) 



Similar results can be obtained for decays of other mesons M + = K + , D + , D^~,B^~ . The coefficients 
B for various decays are given in Table IIII1 The present bounds on the PMNS mixing elements 
\Uni\ for heavy Majorana neutrinos (mjv > 100 GeV) are 23] 



J2 \ Un ^ = ( S l) 2 - °- 005 ' ( S l") 2 - °- 002 ' ( S L ) 2 - °- 010 • 



(19) 



N 



So, Eq. (TT7|) and the bounds in Eq. (fTU|) allow us to interpret the B values in Table IIIII as upper 
bound estimates of the corresponding branching ratios, for the case of heavy Majorana neutrinos 
with masses above 100 GeV,. 

As expect, on r resnhs in »« coincide with those of Re, Q, with discrepancies within 
10% due to variations in the input parameters (a bit larger discrepancies are found in D + and D s 
decays due to the different values we used for fjj and /d s ). The conversion of their theoretical 
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estimates into our notation is a simple factor 10 -14 due to different units used. In the cases where 
Ref. [13] quotes a range, our agreements are with their central values. The sole exception occurs 
in B + — > it~t + t + , where their result is almost exactly a factor 10 larger. We can attribute this 
discrepancy only to a misprint in the power of 10 of their result. 

From Table ITO1 we see that the highest upper bounds (~ 10 -20 ) are for the branching ratios of 
B^ —7- D~£ + £ + , with £ = e or [i. Yet, these bounds are several orders of magnitude too small to 
be detected in current and future experiments, like those at the LHC-6 and Super flavor-factories. 

As a final remark in this subsection, we want to compare these results with what would be 
expected if lepton flavor violation came from other sources, namely Supersymmetry with R-parity 
violation (RPV), or Left-Right symmetric electroweak theories. One can anticipate that such 
LNV processes should involve mass scales well above the electroweak scale, typically around the 
TeV scale, so they can be compared with the same LNV processes mediated by heavy Majorana 
neutrinos. Even though the experimental observation of the LNV meson decays would strongly 
support the hypothesis that neutrinos are Majorana particles, these other sources could produce 
the same signals without involving Majorana neutrinos directly, just as it occurs in neutrinoless 
double beta decays. 

In a supersymmetric extension of the Standard Model that includes RPV, the exchange of 
charged lepton or quark superpartners and neutralinos or gluinos rather than VF-bosons and Majo- 
rana neutrinos can also induce these LNV meson decays. RPV supersymmetry allows for additional 
trilinear terms in the superpotential, of the form: 

W = XijkLiLjEZ + KjkLiQjDl + X'i jk U9Dp c k , (20) 

where i,j,k denote the families, L and Q denote lepton and quark iso-doublet chiral superfields 
and E C ,U C and D c charged lepton and quark iso-singlet chiral superfields. Of these terms, only 
the second leads to LNV meson dec™ ft The effective Logrongian for these decays indnced by 
the RPV terms can be written as [26J 



r 2 

" ' " (21) 



^ 2 = 2^(l + 7 5 )e 



rjPsJpsJps — jVtJt' JT(j,v 



where the hadronic currents are Jps = u a (l + 75)^ and Jj, u = u a a fl ' / (l + 75)^; with color index 

a and a^ u = (i/2)[7 /i ,7^1. Here m v is the proton mass and the explicit forms of the parameters 

P 
r]ps and tjt are given in [27]. We should add that, in the case of LNV decays induced by heavy 

neutrinos, the effective Lagrangian can also be put in the form above if we use r/jy Jy& JvAfi, where 

Jy U A = u a 7 M (l — 75)d Q and t]n = m-l / 10 m GeV ( see sec tion II. c and Appendix C). If we assume 
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that the contribution induced by either gluinos or neutralinos is dominant over the others and the 
masses of the sfermions are almost equal, the parameters rjps and tjt are of the order of 

n hjk^ifk> ( ' o. s m v { a 2 m p \ 

1 

where a s , a 2 , m x , m g and m? denote the strong coupling, SU{2) weak coupling, neutralino mass, 
gluino mass and sfermion mass, respectively, and A^ fe A^ /fc/ actually depends on the process 
(A' 123 A / 111 for B — > Dee, A' 113 A / 112 for B — > Kee, A' 113 A' 111 for B — > -nee, A' 122 A / 111 for D — > Kee, 
A'mAiii fo r D — > nee, A' 112 A' m for K — > -Kee). Besides 0vf3f3 decay, the electron electric dipole 
moment experiments lead to the most stringent bounds on single A' m , which are 5.5 x 10~ 5 for 



Gp-nr- \ 6m g 2m x J 



29l | . There are also several bounds on single 



mf = 100 GeV and 2.4 x 10~ 7 for mi = 1 TeV [28|, 
RPV couplings A' llfc and A' 12fc coming from experimental results for forward-backward asymme- 
tries in the fermion pair production reactions measured at LEP and SLC, and leptonic ir decays, 
respectively; A' llfc < 0.02 and A' 12fc < 0.21 for m rf - = 100 GeV [29]. Imposing those bounds on 
A'ni) A' llfc (fc = 2,3) and A' 12A ,(&; = 1,2,3), the magnitudes of the corresponding terms in Eq. ([22]) 
are of order 10 -8 — 10 -9 for m~ ? = 100 GeV. These should be compared with the parameter 
r]N for the case of LNV induced by Majorana neutrinos. For |£7tv^| 2 ~ 10~ 2 and tun ~ 100 GeV, 
the magnitude of r]]\r becomes unity. Thus, the above bounds on A^ - k would imply that the RPV 
supersymmetric contribution to the corresponding LNV decays should be much smaller than those 
induced by heavy neutrinos for m^ = 100 GeV (and maximally allowed |£7a^|). Conversely, the 
LNV meson decay experiments can be used to put bounds on the corresponding A^ - k parameters, 
especially in those cases where the bounds are very loose or still non-existent. 

Alternatively, a Left-Right Symmetric Model also involves a large mass scale which may char- 
acterize LNV mediated by heavy physics [30j]. In SU{2)i x SU{2)r x U(1)b-l, this gauge group 
breaks down to SU{2)l x U(1)y via an extended Higgs sector containing a bi-doublet <E> and two 
triplets ^l,r whose leptonic couplings generate Majorana neutrino masses and thus lepton num- 
ber violation. The A^^— lepton interactions are not suppressed by lepton masses and have the 
structure L ~ hij Aj"^ Z?(l ± 75)^ + h.c, where the couplings h^ are in general diagonal and 
associated with the heavy neutrino mixing matrix. In this model, short-distance contributions to 
LNV decays arise from the exchange of both heavy right-handed Majorana neutrinos and Al,r, 
which can be parameterized by [311 ] . 



g| 1 gj hjj 

/If 4 M ' M 3 M 2 ' ^ ' 

W R wlv R m W R 1V1 A 

where g 2 , Mw R , Ma and M VR denote the weak gauge coupling, the SU{2)r gauge boson mass, 
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the triplet scalar mass and the right-handed neutrino mass, respectively. These terms are to be 
compared with Tj^Gp/nip corresponding to the LNV decays induced by heavy neutrinos. Imposing 



the current lower bound of 715 GeV on Myy 



n m 



and taking M UR ~ Ma ~ 1 TeV, those terms 



multiplied by m p /G 2 F are of order of 10 9 — 10 10 , which are again very small compared with 
rf$ ax ~ 1 for m N = 100 GeV. 

III. SUMMARY AND CONCLUSIONS 

We have studied lepton number violating decays of charged K, D, D s , B and B c mesons of 
the form M + — > M'~£ + £ + , induced by the existence of Majorana neutrinos. These decays violate 
lepton number by two units, and therefore can occur only if neutrinos are of Majorana type. 
The decays are sensitive to neutrino masses and lepton mixing, and can also provide information 
complementary to neutrinoless double beta decays. We explore neutrino mass ranges tun from 
below 1 eV to several hundred GeV. 

The decay rates are dominated by different weak amplitudes, depending on the mass of the 
neutrinos involved in the intermediate states. 

If the mass of the neutrino that dominates the process is below the mass of the produced 
meson, we find that the main contribution to the branching ratio should come from a two-particle 
intermediate state that goes on shell, formed by a meson (with the correct flavor) and the neutrino. 
These cases have a topology similar to neutrinoless double beta decay. However, the branching 
ratios obtained in these cases are far too small to be detected in foreseen experiments. Indeed, if 
the neutrinos involved are standard (masses below 1 eV, albeit Majorana) we find the branching 
ratios to be below 10~ 31 , and if they are heavier (up to the order of 100 MeV), the branching ratios 
to be below 10 -21 , which are too small to be detected in the foreseen future, so we do not go into 
more refined calculations in these cases. 

Instead, if the neutrino mass is in the range between the masses of the initial and final meson, 
the process is dominated by an intermediate state with just the neutrino, which goes on shell. In 
this "long distance" process, the neutrino is essentially produced and then it decays. Some of the 
branching ratios in this case are now within or near the reach of current or foreseen experiments, 
as shown in Table HTl For example, B(K + — > 7r~e + e + ) can be up to 10 -6 and B(B C — > 7r~e + e + ) up 
to 10 -11 . Experimental exploration of these decays can then at least provide upper bounds for the 
lepton mixing elements of the standard charged leptons with exotic Majorana neutrinos, bounds 
that will be dependent on the mass of the neutrino involved. 
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Finally, if the process is dominated by a neutrino that is considerably heavier than the decaying 
meson, the branching ratio is again far too suppressed to be experimentally observed, as shown in 
Table Hill and previously predicted in Ref. [13]. Indeed, for neutrino masses near 100 GeV or above, 
the branching ratios are all below 10 -20 . In this case we also explore other underlying physics 
sources that could induce these LNV decays without involving Majorana neutrinos directly, 
namely RPV supersymmetric models and left-right symmetric models. Concerning RPV SUSY, 
if we impose the current bounds on the relevant parameters A^ fc , the effect of these interactions 
on the decays would be lower than the effect of neutrinos. Otherwise, in general the experimental 
bound on each of the decays will impose bounds on its corresponding A' parameter. Finally, 
concerning left-right symmetric models, their effect seem to fall far below the contributions of 
heavy Majorana neutrinos, and thus these decays may not be useful to put bounds on those models. 
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Appendix A: The case of Light Neutrinos 

In this appendix we present the calculation of the decay rate for the process M + — > M'~l + £ + 
in the case in which neutrinos are lighter than the mesons in the process. The calculation is done 
assuming that the transition matrix element can be approximated by its absorptive part, which 
given in Eq. Q. 

After doing the integral in Eq.(j4|), the square of the amplitude, summed over the final lepton 
spins, becomes: 

\M abs \ 2 = ^\VcbV ud \ 2 F+ 2 D F+ 2 D \U 2 m rn N \ 2 ^- x T, (Al) 

where we have defined T = 

2Tr [l/ 2 {i> D i B + m 2 D0 + £ D o( 7 ° J*b+ ^7°)) ft {i>B i>D + m 2 D0 + £ D o(^7° + 7° M)] , 

and where |pjv| is the neutrino 3- momentum in the rest frame of the D-N pair. Using the well 
known expression 

\(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2yz - 2xz, (A2) 

it can be written as |pjy| = X 1 ' 2 (rn 2 De ,rn 2 ~ )o ,'m 2 N )/2rnDe- 

Since the expression for|A^ a b s | 2 in Eq. (|A1|) is not explicitly covariant, it is convenient to separate 
the phase space integral over the D-£-£ final state (dps 3 ) into the 2-body phase spaces for B — > 
t\ + Xd£ and Xf)£ — > D + £2, with the invariant mass of the pair Xd£ integrated over its physical 
range: 

/I 3 S'v- f dm 2 !' f 

dps 3 = J H ^-^_(2vr) 4 5 4 (S^ -p M ) = J -^ J dps (B ^ lXm) J dj>s {Xm ^ Dh) , (A3) 

where the 2-body phase spaces in their respective rest frames reduce to: 

dpS {B ^ lXm) = Y^2^ dQ ^ d P S (X D ^Dl 2 ) = Y^2^ dn ^ 

and the 3-momenta in the respective cases are: 



~ X^ 2 (m 2 Bi m 2 Dt ,m 2 ) ^ \V 2 (m 2 DV m 2 D ,m 2 ) ^ 

2niB 2rri£)e 

Now, the integration over df^ 2 of the non trivial factor T in Eq. (|Aip can be expressed as: 



(A4) 



/ 



dVS {Xm ->B h) T = j^^L 4nU, (A5) 
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where 1Z is a long expression of dimension m . 



K = | 8(m^ + 2E z D0 yE l E 2 + l6(m z D o + 2E Z D0 )(E D E 2 + \\*\*){E B Ex - \h\l (A6) 

+ W(m 2 D o + 2E 2 DO )E D o ( E 2 (E B E 1 - \h\ 2 ) - E 1 (E D E 2 + |1 2 | 2 ) 



+ 16E D o ( 2E B {E B E l - \\ 1 \ 2 ){E D E 2 + |1 2 | 2 ) - m 2 D E 2 (E B E 1 - \h\ 2 ) - m 2 B E l {E D E 2 + \l 2 \ 2 ) 

2 / Z7' ZT'/TZT' IT 1 l2 i 1 1 I 2 I || I 2~\ I c\{ rp rp \ / IT 1 1 1 I 2 ZT'll 1 2 \ oil I 2 1 1 |2 



+ 8E 2 D0 { ExEzUEb - E D Y + \h\ 2 + \h\ 2 ) + 2(E B - E D )(E 1 \l 2 \ 2 - E 2 \h\ 2 ) - 2\h\% 
+ 8(m 2 D E 2 - 2E 2 D E 2 - 2E D \\ 2 \ 2 ){m 2 B E l - 2E 2 B E 1 + 2E B \h\ 2 ) }, 

where all kinematical variables here are defined in the rest frame of the D-l pair and are functions 
of its invariant mass mjjf. 

Ed = m m + m 2 D -m 2 ^ = m 2 m - m 2 D + m 2 = m 2 m + m 2 D0 - m 2 N 

2m m ' 2 2m De ' D 2m Di 

m\ + m 2 m - mj m| - m 2 m - m 2 A^m^m^mf) 

Hi B — , Hj\ — , | i^i — . 

2m D £ 2m B i 2m m 

Finally, since the result in Eq. (|A5P is independent of angles, the subsequent integration over d£li x 
simply brings a factor Air. The decay rate T(B + — > D^l^i^) then results in the expression: 

(rriB-rrn) 
tVp+ , T^-n+o-\_ G F it/ ^ |2i7+2 F +2 \U 2 Nt m N? f dm Di \p N \ 2 - 

T{B -+DU)- j^IVM F BD F DD ^ J ^— -Hm |1 2 | xK, (A8) 

(mjj+mi) 

where |pat|, |1i|, II2I an d 1Z were defined above and are explicit functions of m B £. The integral in 
the expression above can be easily done numerically. 



Appendix B: The case of intermediate mass neutrinos 

Here we present the meson decay rate M + — > M'~£ + l + in the case in which neutrinos have a 
mass in the intermediate range tum 1 < ttin < rriM- 

The square of Ai in Eq. ([9]) and sum over external lepton spins of this reduced amplitude, after 
some algebra, can be written as : 

\M\ 2 = 32m 2 N \{m 2 N - m]) 2 {h ■ l 2 ) + mj ((m 2 N - mj) 2 - m 2 M m 2 M ,) } (Bl) 

The final 3-body phase space can again be separated into two 2-body integrals, and another 
over the invariant mass of the intermediate state (which in this case is the neutrino momentum 
squared, p%)' 

/ dps 3 = / -^- I dps {M ^ hN) / dps {N ^ hM ,), (B2) 
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where dpS( M ^ h N ) = 1/(167t 2 )(|1i|/tom)^i and dpS( N _+ h M ,^ = l/(167r 2 )(|l2|/mAr)cftl 2 . This time, 
the propagator of the intermediate neutrino in the matrix element [see Eq. ([8]) ] can be approximated 
by a delta function, since it is a narrow state: 

^ K( 2 2 \ 

T~2 9 \9 , 2~n2 ^ F "VR/V ~~ m AfJ- 

iPtf - m N) + m ^ r Ar "^IV 

The only term in \ftA\ 2 that depends on an integration angle is the one that contains {l\ ■ I2) [see 
Eq. (jBip ]. In the neutrino rest frame, J dQ.2{W • h) = ^E\E2, where E\ = {m 2 M — m 2 N — mf)/2m^i 
and E2 = {m% + rn 2 — m 2 M ,)/2mN are the respective energies of the external leptons. All other 
solid angle integrals give just a factor 4ir. Putting everything together, we obtain the decay rate: 

F(M -> MW) = ^^!l t fl A V, Q V n ^\U m \* J i VH ^—^—^ (B3) 

Gf f 2 f 2 it, T/ |2 l^iwl |ll| l^l 



-/m/m'I^jQ^; 



327r 2 m M J M ' M ' ' '" V 92 91 ' m N T N m M m N 

x {{m 2 N + mf)m 2 M - (m 2 N - mj) 2 } {{m 2 N - m)) 2 - {m 2 N + mj)m 2 M ,} , 

where |li| = \ 1 ' 2 (m 2 M ,m 2 v ,rn 2 )/2rnM and | X^ I = A 1 ' 2 (m 2 v ,m 2 v/ /,m,|)/2m,Ar. 



Appendix C: The case of heavy neutrinos 

Here we present the derivation of the meson decay rate for M + — > M'~l + £ + (£ = e, /u, r) when 
the exchanged Majorana neutrino is heavier than the decaying meson. In this case, both weak 
amplitudes shown in Fig. 1 contribute with similar strength to the process. 

In order to see the approximations involved, let us first consider the neutrino mass m^ not to 
be necessarily high. Let us first recall the valence quark content of the decaying meson M as (Qq) 
and of the produced meson M' as (9291)- Let us also denote by J,^ . = 07^(1 — 75)5 the weak 
V — A quark current with flavor change q —> Q. In much the same way as before, we rearrange 
the lepton line using charged-conjugated spinors and the Majorana character of the neutrino, thus 
appearing an irrelevant Majorana phase Atv = exp^Ar). Then the contribution of Fig. l.a (and 
crossed diagram) can be written as: 

Mu = (-1) ^ u&Xi, (v q * lQ v; q2 ) (M'\j^ 2q) r {Qqi) \M) (ci) 

-n\ r-i \( ¥N+m N ¥N + m N \/i \ n\ 

x ujihh^l + 75) 2 m 2 T-fT— - + T7-2 2—- 7,(1 " 75)^(^i)- 

y K N - m N + 11 NrriN k N - m z N + iT N m N J 



Here we called k^ and k N the corresponding neutrino momenta for the two crossings of the external 



lepton lines. Now, if the neutrino is heavy, in the denominators we neglect all except m 2 N . In 



23 



addition, using the orthogonality of the chiral projectors and the relation {7^,7^} = 2g^ u , the 
expression reduces to: 

Mia = ^f U$ t X *n K Q V* q2 ) (M'\J^ q) J {Qqi)ll \M) ^{ui(h){l - TBMh)]. (C2) 

In much the same way one can treat the contribution of Fig. l.b, but now the roles of the flavors 
q and q\ are interchanged, thus the quark currents and CKM elements are different: 



M 



n 2 

^F TT*2 X* 



u = -f UffiX*N {V q * Q V q * iq2 ) (M'\J^ qi) J {Qq) ,\M) —[ui(h)(l - -ysMh)]. (C3) 

The largest uncertainty here is in the determination of the hadronic matrix element. As a first 
estimate, reasonable for the level of accuracy we seek, we can separate the currents inserting the 
vacuum and assuming it saturates the expression (vacuum saturation approximation). We must 
then do a Fierz rearrangement in M\ a to match the flavors, which then mismatches the color, 
inducing a suppression factor 1/N C . The hadronic currents then reduce to their corresponding 
decay constants, /m and /m', and the total amplitude for heavy neutrino exchange ("h") becomes: 



My 



Mla+M lb 

G 2 

2 



r F TT*2 \* 

u Nl A N 



v*„v* 

T /* T /* 1 giQ <?<?2 



qQ qiQ2 



N r 



fMfM'iPM ■ PM>) M^)(l - 75)^(Zl)]. ( C4 ) 

m N 



One may neglect the 1/N C term, except if the other term is much more CKM-suppressed. The 
square and sum over final polarizations of this amplitude is: 



\M h \ 2 = |/C h | 2 32 {pM-PM'f (h-h), 



(C5) 



where we have gathered all constant factors under the symbol 



l/Ci 



1 2 _ n^ 
— t*F 



TT*2 



m N 



V* „V* 
V* V* 4- qi ® qq ' 2 



ImIm 1 - 



The decay width T(M + — > M' £ + £ + ) can now be calculated explicitly, by integrating |A^h| 2 ° ver 
the final phase space. This time we express the M'-t-t phase space as: 

dm?no 



dps 3 



2tt 



^P S (M^M' x u ) \ dvs(x u -+u) 



where dps (M ^ M/ Xu) = (1/167t 2 )(|p M 'IMm)^m' and dps^ Xu ^ u) = (l/l6ir 2 )(\l 2 \/m ee )dn £2 . In 
the frame of the lepton pair, the integrand (|C5P is independent of angles, so the integral over d£le 2 
is a simple factor of 4ir, and equally for the subsequent integral over d&lM'- 
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T(M+ -* M"tt) = IV^ f ^^^(rnl + mlt,- m\f {ml - 2mj), 
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The decay rate is then: 

(m M -m M ,) 2 

2! 8tt 2 itim J 27r tum ram 

(C6) 
where \pm'\ = A 1 ' 2 (m^ f ,m^ f/ ,r7T,^)/2?TT,M and 1 1^ I = A 1 ' 2 (m^,m|,rri 2 )/2ma. The factor 1/2! 
appears because there are two identical particles in the final state. This integral can be easily done 
numerically. 
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